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摘要 :针对 干旱 牧区 生态 脆弱 性 特点 ,以 GIS 为 平台 ,利用 气象 .土壤 、 土 地利 用、 水 资源 开发 利用 等 数据 ,结合 SWAT 模型 模拟 
结果 ,建立 了 一 种 基于 SWAT 模型 子 流域 划分 的 干旱 牧区 生态 脆弱 性 评价 体系 , 对 西北 干 量 牧 区 艾 布 盖 河 流域 1990 年 .2010 
年 生态 脆弱 性 进行 了 评价 。 从 评价 结果 可 以 看 出 ,研究 区 北部 及 西部 生态 稳定 性 保持 较 好 ,生态 脆弱 等 级 下 降 主 要 发 生 在 研究 
区 的 东南 部 ,20 年 间 耕 地 扩张 造成 的 草地 被 侵占 及 水 资源 过 度 开 发 利用 是 造 威 华 态 退化 的 主要 原因 。 评 价 体系 紧 紧 围绕 干旱 
牧区 “水 - 草 - 畜 ” 关 系 设 置 评价 指标 ,将 观测 数据 与 水 文 模型 模拟 相 结合 ,统计 数据 与 空间 地 类 相 联系 ,点 面 结 合 ,能 够 较 好 的 反 
映 出 研究 区 20 年 间 生 态 脆 弱 程度 变化 趋势 及 主要 问题 。 
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Abstract: Ecological vulirerability evaluation is a prerequisite to conducting ecological restoration. The common ecological 
vulnerability evaluation methods have some applicability problems in arid pasturing areas: 1) the evaluation indicators are 
numerous, and the datalare Usually difficult to obtain; 2) an administrative region, such as a city or a country, is often 
used as the evaluation unit, which is too large and lacks specific significance. Commonly used data, such as meteorological 
data, soil data, land-use data, and water resource development and utilization data, as well as the simulation results in the 
SWAT niodel, are used in this research to build an ecological vulnerability evaluation of arid pasturing areas based on the 
SWAT model Ecological vulnerability of arid pasturing areas has been divided into ecological vulnerability of the natural 
ecosystemMand ecological vulnerability of the human-nature coupling system. Drought index, vegetation coverage, soil bulk 
density ,Jand evapotranspiration capacity of soil water are used to depict ecological vulnerability of natural ecosystem ; 
evenness index, contagion index, ecological suitability index, and water utilization modulus are used to depict the human- 
nature coupling system. Firstly, the SWAT model has been used to divide sub-basins and hydrological simulations by the 


combination of observed data and simulated data, and the association of the statistic data with spatial land type. All 


基金 项 目 : 国 家 自然 科学 基金 项 目 (51579157) ; 中央 分 成 水 资源 费 项 目 (1261430122054) ;中 国 水 利水 电 科 学 院 科研 专项 项 目 ( MK2016J16) 
收 稿 日 期 :2016-04- 19; 网 络 出 版 日 期 :2017-02-22 
* 通讯 作者 Corresponding author. 上 -mail : songyf_90@ 163.com 


http :// www.ecologica.cn 


3806 生 态 学 报 37 卷 


indicators have been set based on CIS linear piecewise membership functions have been adopted to calculate index values, 
and analytic hierarchy process (AHP ) is adopted to determine the weighs of criteria layers and index layer. The ecological 
vulnerabilities of the research area are divided into five levels. This evaluation model was used in the Aibugai Basin (which 
is located in arid pasturing areas in Northwest China) to evaluate the change in ecological vulnerability caused by land-use 
variation between 1990 and 2010, and 17 sub-basins were divided according to demand. The results showed that, during a 
period of 20 years, both the entire research area and each sub-basin went through a CIEV (comprehensive index lof 
ecological vulnerability) decrease to some degree. Sub-basin 9, 10, 11, 12, and 15 developed a moderate vulnerability 
level from a slight vulnerability level in 1990, and sub-basin 16 and17 have changed to a sever vulnerability level from a 
moderate vulnerability level in 1990. The CIEV of the entire research area decreased to 0.59 in 2010 from 0.70 in 1990. The 
northwestern part of the research area had better ecological stability, and the degradation in ecological Vulnerability levels 
took place mainly in the agro—pastoral ecotone located in the southeast of the research area. With thé expanding demands of 
grain and cattle, agricultural areas have gradually expanded northward in the past 20 years. Declines in grassland areas and 
overdeveloped water resources caused by farmland expansion were the main reasons that causéd the eeological degradation in 
the research area. Closely tied to the relationship with water, grass and cattle in arid pasturing/areas, to establish indicators, 
this evaluation model effectively reflected the aggravation of ecological vulnerability cauSed by czopland expansion in the 


research area from 1990 to 2010. 
Key Words: arid pasturing areas; SWAT; ecological vulnerability; evaluation method 


生态 脆弱 性 评价 是 了 解 区 域 或 流域 生态 脆弱 问题 与 合理 弄 展 生态 修复 的 前 提 '"”。 我 国 脆弱 性 研究 开 
台 于 20 世纪 90 年 代 “ ,研究 区 域 由 最 早 的 西南 喀斯特 典型 脆弱 区 “逐渐 扩展 到 北方 农 牧 交错 带 “! 南方 丘 
陵 区 '"| 青藏 高 原 区 “” , 近 些 年 ,西北 干旱 . 半 干 时 牧区 的 生态 脆弱 性 问题 逐渐 为 人 们 所 关注 “1 。 然 而 ,目前 
常用 的 脆弱 性 研究 方法 在 我 国 西北 干旱 牧区 仍 唱 存在 着 一 些 适 用 性 问题 ,我 国 西北 干旱 牧区 植被 类 型 单一 ， 
主要 为 草地 ,社会 经 济 结构 相对 简单 ,以 草地 农 牧 业 为 主 ,工业 化 程度 较 低 , 这 种 情况 下 ,过 于 庞杂 的 指标 体系 
增加 了 数据 获取 与 指标 计算 难度 的 同时 s 也 在 一 定 程度 上 降低 了 评价 体系 的 针对 性 ,此 外 ,西北 干旱 牧区 地 ) 
人 稀 , 各 类 测 站 和 统计 资料 支持 能 力 较 为 有 限 , 较 常用 的 行政 区 域 评价 单元 往往 范围 过 大 而 缺乏 局 部 指导 
综 上 所 述 ,本 文 以 西北 毛 早 牧区 艾 布 关 河 流域 为 例 ,以 GIS 为 平台 ,利用 气象 土壤、 土地 利用 水 资源 开 
发 利用 等 数据 资料 ,结合 SWAT 分 布 式 水 文 模型 ,建立 了 一 种 基于 SWAT 模型 子 流域 划分 的 干旱 牧区 生态 脆 
弱 性 评价 体系 。 


1 评价 体系 构建 方法 


本 交 将 干 星 牧区 生态 脆弱 性 划分 为 自然 生态 系统 脆弱 性 及 人 类 -自然 耦合 系统 脆弱 性 ,其 中 自然 生态 系 
统 反映 盯 自然 生态 本 底 或 基底 ;人 类 -自然 耦合 系统 脆弱 性 反映 了 人 类 社会 系统 与 自然 系统 间 的 相互 作用 ,是 
人 类 活动 对 自然 生态 系统 产生 影响 的 体现 。 
1. 和 基于 SWAT 模型 的 评价 单元 划分 

SWAT(Soil and Water Assessment Tool) 是 一 款 具 有 很 强 物理 机 制 的 分 布 式 水 文 模型 ,能够 模拟 多 种 土壤 、 
土地 利用 和 管理 条 件 下 的 水 文 过 程 (521 。SWAT 模型 由 SWRRB 模型 演化 而 来 ,吸收 了 SWRRB 模型 划分 子 
流域 的 功能 ,能 够 基于 DEM 数据 进行 子 流域 划分 ,本 文采 用 SWAT 模型 对 干旱 牧区 艾 布 盖 河 流域 进行 水 文 模 
拟 , 并 在 水 文 模拟 的 基础 上 以 子 流域 为 评价 单元 建立 了 干旱 牧区 生态 脆弱 性 评价 体系 。 这 种 基于 分 布 式 水 文 
模型 子 流域 划分 的 评价 方法 在 我 国 西北 地 区 更 具 灵 活性 与 适用 性 ,主要 优点 如 下 :(1) 每 个 子 流 域 作为 独立 
的 产 汇流 单元 具有 近似 的 水 文 特征 ,在 水 文 计算 中 更 具 合理 性 ,便于 进行 流域 综合 规划 管理 ; (2) 子 流域 基于 
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DEM 数据 生成 ,可 根据 需要 在 SWAT 模型 中 通过 设置 子 流域 临界 集 水 面积 CSA( Critical Source Area ) 调整 其 
大 小 ,更 具 灵 活性 ; (3) 在 分 布 式 的 SWAT 模型 中 ,每 个 子 流域 单独 进行 水 文 过 程 计算 ,模拟 结果 可 用 于 评价 
指标 计算 ,以 弥补 西北 干旱 地 区 数据 资料 匮乏 的 问题 。 

1.2 评价 指标 选取 与 计算 方法 

选取 恰当 的 评价 指标 对 “水 - 草 - 畜 ”关系 进行 合理 刻画 是 干旱 牧区 生态 脆弱 性 评价 的 关键 ,本 评价 体系 选 
取 能 够 反映 气候 特征 的 干旱 指数 ,能 够 反映 植被 条 件 的 植被 盖 度 ,能 够 反映 土壤 属性 的 土壤 容重 ,以 及 能 够 反 
上 映 天 然 植被 生态 需 水 情况 的 土壤 水 薰 散 能 力 这 4 个 指标 对 研究 区 自然 生态 系统 脆弱 性 进行 刻画 。 

人 类 活动 通过 改变 下 垫 面 对 生态 系统 施 以 影响 ,土地 利用 变化 是 人 类 活动 与 自然 生态 系统 相互 代用 最 直 
观 的 体现 ,土地 利用 变化 可 以 通过 景观 格局 加 以 反映 。 此 外 ,每 种 土地 利用 类 型 都 承载 了 一 定 的 生态 功能 价 
值 和 水 量 需求 ,使 其 在 一 个 地 区 长 久 以 来 形成 的 自然 环境 中 具有 特定 的 适宜 性 和 水 资源 脆弱 性 4 本 文选 取 能 
够 反映 景观 格局 变化 的 均匀 度 . 草 延 度 ,能 够 反映 土地 利用 类 型 生态 价值 功能 和 稳定 性 的 生态 适宜 指数 ,以 及 
反映 水 资源 开发 利用 水 平 的 水 资源 利用 模 数 4 个 指标 对 人 类 -自然 耦合 系统 脆弱 性 进 簿 刻画 图 1) 。 
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图 1 干旱 牧区 生态 脆弱 性 评价 体系 


Fig.1 A ‘The evaluation system on ecological vulnerability of arid pasturing areas 


1.2.1 自然 生态 系统 脆弱 性 指标 
干旱 少雨 的 气候 条 件 是 导致 我 国 西北 干旱 牧区 生态 脆弱 的 首要 原因 ,干旱 牧区 降水 具有 极 强 的 时 空 变异 
性 ,在 强烈 的 蒸发 作用 下 ,降水 的 分 布 很 大 程度 上 决定 了 生态 系统 的 脆弱 性 程度 ,因此 ,本 评价 体系 采用 干旱 
此 数 D1 反映 各 予 流域 当年 的 干旱 程度 : 


D1 = PET, + PRECIP， (1) 

式 中 每 个 洱 流 域 的 潜在 蒸 散 量 PET, 和 降水 量 PRECIP, 可 在 SWAT 模型 计算 结果 中 直接 读 出 。 

我 国 西北 干旱 牧区 地 带 性 植被 多 为 旱 生 丛生 草本 植物 ,植被 盖 度 常 被 用 来 进行 草地 退化 程度 判别 ,是 干 
量 地 区 植被 生长 及 分 布 情况 的 直观 体现 。 本 评价 体系 植被 盖 度 指标 计算 先 由 遥感 数据 反 演 得 到 高 中 低 盖 度 
草地 类 型 ,再 对 每 个 子 流 域 进行 面积 加 权 得 到 。 

土壤 尤其 是 根系 层 土壤 作为 “水 - 草 ” 相 互 作用 的 重要 介质 ,是 生物 , 非 生 物 间 物质 能 量 交换 最 为 活跃 的 
场所 ,其 属性 优 劣 也 是 生态 环境 的 综合 反映 己 ] 。 土 壤 容重 与 土壤 质地 、 级 配 土壤 颗粒 密度 .土壤 有 机 质 含 
量 、 和 孔隙 度 等 因素 有 关 , 是 土壤 重要 的 物理 参数 之 一 。 一 般 来 讲 ,容重 较 小 的 土壤 ,结构 政 松 ,级 配 良好 ,有 机 
质 含量 高 ,有 利于 土壤 通 水 通气 。 由 于 干旱 牧区 草本 植被 根系 埋 深 较 浅 ,多 分 布 在 30 cm 以 内 ,本 评价 体系 中 
土壤 容重 指标 采用 FAO 土壤 数据 中 0 一 30 cm 土壤 层 属性 数据 在 每 个 子 流域 面积 加 权 得 到 。 
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根系 层 土壤 水 是 天 然 植被 生态 需 水 的 直接 来 源 , 同 时 也 为 其 返青 阶段 提供 了 必要 的 水 源 储备 ,在 天 然 植 
被 生长 发 育 的 各 个 过 程 中 起 着 至 关 重 要 的 作用 。 在 干旱 牧区 ,大 部 分 降水 以 土壤 薰 发 和 植被 莱 腾 的 形式 
消耗 (实际 蒸 散 量 £7) ,用 以 维系 地 带 性 植被 生长 '* ,本 评价 体系 采用 土壤 水 蒸 散发 能 力 指标 SETC, 反映 各 
子 流域 土壤 水 分 可 维系 天 然 植被 蒸 散 发 的 能 

SETC, = SW, + ET, (2) 

式 中 ,各 子 流域 土壤 水 量 SW, 和 实际 蒸 散 量 57, 可 从 SWAT 模型 计算 结果 中 直接 读 出 。 
1.2.2 ”人 类 -自然 耦合 系统 脆弱 性 指标 

均匀 度 ( SWE, ) 是 反映 各 种 景观 类 型 相对 面积 比重 的 指标 5 , 取 值 范围 为 0 一 1, 若 各 种 景观 类 型 面积 比 
例 差异 越 大 ,其 值 越 趋 近 于 0; 各 类 型 面积 比例 越 接 近 , 越 趋 近 于 1。 其 计算 公式 为 : 


= 2p, Xx lnp， 
SWE, = A (3) 


lInm 


__，。 式 中 ,为 第 j 种 类 型 的 景观 在 整个 景观 中 所 占 的 面积 比例 ;m 为 景观 类 型 总 数 允 于 我 图 西北 干旱 牧区 景 

> ” 观 类 型 以 草地 为 主 ,其 他 景观 类 型 面积 远 小 于 草地 面积 , 则 可 认为 均匀 度 越 小 , 疆 地 爱 风 类 活动 影响 越 小 , 反 

0 入 ”之 越 大 。 该 指标 由 每 个 子 流域 土地 利用 数据 计算 得 出 。 

< 草 延 度 (CONTACG ) 是 反映 景观 斑 块 联通 程度 和 延展 趋势 的 一 个 指标 , 取代 为 0 一 100, 营 延 度 越 高 ,说 明 

景观 中 的 某 一 种 优势 斑 块 类 型 具有 较 好 的 连通 性 ;反之 ,说 明 景观 臣 砚 化 程度 较 高 ,景观 格局 继续 变化 的 风险 

,一 ”大 。 这 一 指标 通过 顶 格 分 析 计算 ,利用 各 子 流域 的 土地 利用 栅 格 图 在 .PRACSTATS 软件 中 计算 得 出 。 

©O 地 带 性 景观 是 一 个 地 区 漫长 演变 过 程 中 逐渐 形成 的 最 适 状 态 ) 生态 适宜 指数 表征 了 不 同 景观 类 型 在 当地 

| 一 生态 系统 中 的 适宜 程度 ,该 指数 基于 不 同 土地 利用 类 者 生态 芒 能 及 适宜 程度 ,通过 打分 给 出 。 该 指数 取 值 

@D ”0 一 100, 值 越 高 说 明 该 种 景观 类 型 在 当地 的 生态 系统 中 芍 活 宜 性 越 强 ,需要 说 明 的 是 我 国 西北 干旱 牧区 耕地 

| 中 包含 大 量 灌溉 饲 草料 地 ,能 够 为 牲畜 提供 饲料 来 源 , 春 一 定 程度 上 有 利于 草地 的 涵养 ,因此 生态 适宜 性 较 高 

于 居 工 地 ;林地 多 为 非 地 带 性 植被 ,因而 略 低 于 草地 \ 具 体 打分 见 表 1。 指 标 值 计算 根据 不 同 土地 利用 类 型 生 
态 适 宜 性 指数 在 每 个 子 流域 面积 加 权 得 到 。 


表 1 不 同 土地 利用 类 型 生态 适宜 性 指数 表 
Table 1 Index of ecological suitability in each type of landuse 


Toil 


Dd 土地 利用 类 型 耕地 林地 草地 水 域 居 工 地 未 利用 地 
Land use type Farmland Forest Grassland Waters Settlements Unutilized land 
生态 适宜 性 


Ecological suitability 一 中 对 光 

社会 经 济 发 展 水 平 超过 当地 水 资源 承载 能 力 造成 的 “水 - 草 - 畜 ”失衡 是 导致 干旱 牧区 生态 脆弱 的 重要 原 
因 ,因此 可 将 水 资源 开发 利用 情况 作为 干旱 牧区 生态 脆弱 性 程度 的 表征 之 一 ,本 文采 用 水 资源 利用 模 数 
(MWRU,) 来 反映 每 个 子 流域 水 资源 开发 利用 水 平 。 计 算 时 先 算出 每 种 土地 利用 类 型 的 水 资源 利用 模 数 , 即 
用 当年 灌溉 水 量 除 以 耕地 面积 得 到 耕地 类 型 的 水 资源 利用 模 数 、 用 当年 二 三 产 以 及 生活 用 水 量 之 和 除 以 居 工 
地 面积 得 到 居 工 地 类 型 的 水 资源 利用 模 数 ,并 规定 天 然 状态 下 的 土地 利用 类 型 水 资源 利用 模 数 为 0, 最 终 由 
每 种 土地 利用 类 型 水 资源 利用 模 数 在 子 流域 面积 加 权 得 到 各 个 子 流域 当年 的 水 资源 利用 模 数 。 

至 此 ,本 评价 体系 中 所 有 指标 均 建 立 完毕 , 见 表 2。 
1.3 评价 指标 生态 脆弱 性 指数 计算 方法 

为 了 使 每 个 取 不 同 单位 、 量 级 的 指标 值 都 能 介 于 [0,1] 区 间 , 定 义 生 态 脆 弱 性 指数 ( IEV, ) , 且 规 定 生态 
脆弱 指数 越 接近 于 0, 则 该 项 评价 指标 所 代表 的 内 容 生态 脆弱 性 越 差 , 越 接近 于 1 则 越 好 。 评 价 指标 生态 脆 
弱 性 指数 的 计算 采用 线性 分 段 隶 属 度 函 数 法 '"" ,该 方法 将 指标 分 为 正 向 指标 及 逆向 指标 , 正 向 指标 的 正 V， 
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随 指标 值 的 变 大 而 变 大 (图 2a 式 (4) ) ,逆向 指标 的 正 V, 随 指标 值 的 减 小 而 变 大 (图 2b、 式 (5)) ,指标 节 ， 


aocisdie 分 别 对 应 IEV, 数值 的 0.2.0.4.0.6.0.8 1.0。 


准则 层 Criterion layer 


表 2 评价 指标 与 数据 来 源 表 


Table 2 The indexes and data 


指标 层 Index layer 


数据 来 源 Data sources 


点 值 


然 生态 系统 脆弱 性 


Vulnerability of the natural ecosystem 


人 类 -自然 耦合 系统 脆弱 性 


Vulnerability of the human-nature 


coupling system 


生态 脆弱 性 指数 IEV; 
Index of ecological vulnerability 


Ci 


干旱 指数 
植被 盖 度 /% 

土壤 容重 / (g/cm’ ) 
土壤 水 蒸 散 能 


均匀 度 


蔓延 度 


生态 适宜 性 


水 资源 利用 模 数 /( 万 m3/km?) 


1 
1 
1 
1 
| 
1 
1 
1 
1 
d; 


IEVisg = 


Ci 
指标 值 
Index value 
2 IEV; 曲 线 
Fig.2 The IEV; curve 


SWAT 计算 结果 
遥感 数据 反 演 
FAO 土壤 数据 
SWAT 计算 结果 

土地 利用 数据 

土地 利用 数据 

专家 打分 

水 资源 公报 等 水 资源 开发 利 


Xi 
2 = Xi a, 
Ci 
xX.— a. 
0.2 + 0.2| 一 ] a; < x; bb, 
b= 
.= 
0.4 + 0.2 | b, < xi 到 ci 
人 
WV 
0.6 + 0.2 ] ci <%,d, 
太一 ci 
Xi —d; 
0.8 + 0.2 | d, < x; < e; 
0 
1 e, <x 
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1 Xi <e, 
ee 
0.8 + 0.2| 一 el <x,d, 
本 二 已 
= 
0.6 + 0.2 d, < %, Ee, 
d, = 
IEV,ys, = ;一 %; 习 
a 地 六 有 让 & 
多 =, 
b; =w, 
0.2 + 0.2 b, < wx; qi 
0 加 Ci 
Xi 
由 2 一 Xi > Qi 
Ci 


1.4 权重 的 确定 及 生态 脆弱 性 等 级 的 划分 
准则 层 和 指标 层 权 重 的 确定 采用 层次 分 析 法 (AHP) ,各 指标 生态 脆弱 性 指数 IEV, 乘 以 权重 得 到 各 子 流 
域 生态 脆弱 性 综合 指数 ( CIEV, ,Comprehensive Index of Ecological Vulnerailit : 
CIEV,， = > ps IEV, (6) 


式 中 , w，、w 分 别 为 指标 层 和 准则 层 权重 , IEV, 为 第 个 指标 续 态 脆 妮 性 指数 。 
本 评价 体系 将 研究 区 生态 胸 弱 性 划分 为 不 脆弱 . 轻 度 脆 避 中 庆 脆 弱 .重度 脆弱 .极端 脆弱 5 个 等 级 ,脆弱 
性 程度 依次 增强 ,不 同 脆弱 等 级 描述 见 表 3。 


表 3 研究 区 生态 脆弱 性 等 级 划分 表 


Table 3 Division of ecological vulnerability level of the research area 


脆弱 等 级 CIEWE 取 值 范 围 特征 描述 

Vulnerability grade Valire range of CIEVi Characterizations 

不 脆弱 None vulnerability 0.8 三 1.0 生态 环境 稳定 ,景观 格局 合理 ,人 类 干扰 程度 极 低 

轻 度 脆弱 Slight vulnerability 0.6 一 0.8 生态 环境 较 稳 定 ,景观 格局 合理 ,存在 轻微 人 类 活动 干扰 
中 度 脆 弱 Moderate vulnerability 0.4 一 0.6 生态 环境 不 稳定 ,人 类 活动 对 生态 环境 存在 明显 干扰 
重度 脆弱 Severe vulnerability 0.2 一 0.4 生态 环境 受到 人 类 活动 的 严重 干扰 ,景观 格局 不 合理 
极端 脆弱 Extreme vulnerability 0 一 0.2 各 种 生态 环境 问题 突出 ,生态 环境 极端 脆弱 


CIEV,: 生态 脆弱 性 综合 指数 Comprehensive index of ecological vulnerability 


2 实例 应 用 


2.1 研究 因 概 况 

艾 布 盖 河 沉 域 ( 百灵 庙 水 文 站 上 游 ) 主要 位 于 内 蒙古 自治 区 包头 市 达尔 罕 成 明 安 联合 旗 境 内 ,流域 面积 
5192 Ra2( 图 3) 。) 达 尔 罕 茂 明 安 联合 族 北 接 蒙 古国 , 南 邻 呼和浩特 市 的 武川 县 .包头 市 的 固 阳 县 , 西 与 巴 彦 淖 
尔 市 的 乌拉 特 中 族 相 邻 , 东 与 乌 兰 察 布 市 的 四 子 王 旗 交 界 ,是 典型 的 西北 干旱 牧区 畜牧 业 旗 县 。 研 究 区 冬季 
寒冷 千 燥 ,春季 干旱 少雨 ,多 年 平均 降水 量 253.45 mm ,多 年 平均 水 面 蒸发 量 为 2480.57 mm。20 世纪 90 年 代 
以 来 , 随 着 经 济 社会 的 不 断 发 展 ,研究 区 南部 农 牧 交错 带 农田 面积 不 断 扩张 ,景观 破碎 化 程度 进一步 加 剧 , 草 
原 牧区 生态 环境 脆弱 问题 日 益 凸 显 ( 图 4)。 
2.2 SWAT 模型 水 文 模拟 

使 用 基于 ArcGIS 9.3 的 SWAT 2009 版 本 进行 模型 的 构建 。SWAT 模型 输入 数据 见 表 4, 气 象 水 文 站 分 布 
情况 见 图 $ ,空间 数据 统一 使 用 北京 54 坐标 系 。 通 过 设 定子 流域 临界 集 水 面积 CSA 对 子 流 域 进行 划分 "5 ， 
本 例 结合 需要 将 CSA 设 定 为 20000 hm? ,所 得 17 个 子 流域 如 图 6 所 示 。 
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® 苏 木 、 乡 、 镇 
一 ”县 界 
一 - 苏 木 乡 、 镇 界线 
-一 河流 
比例 尺 1:3000000 
109° 110° lIE 
图 3 研究 区 地 理 位 置 
Fig.3 Geographic location of the research area 
800 
< = 1990 年 
过 ”600 2010 年 
党 400 
岩 200 
0 
~ 4000 
点 3000 上 
车 2000 上 
号 1000 四 
哄 0 一 一 一 a mm 
袜 号 竺 多 究 号 餐 58 拒 吕 袜 号 
蓉 呈 长 5 机。 羡 六 虹 5 里 
图 例 景观 类 型 变化 = 四 甫 3 爱 吕 
口 流域 界线 口 无 变化 EE 5 员 关 三 
z 号 草地 转业 地 E 
+ [S] 
110° 111°E 
图 4 景观 变化 情况 图 
Fig.4 Variation of landscape in the research area 
表 4 SWAT 模型 所 需 的 数据 
Table 4 Data used in SWAT model 
数据 名 称 Data 数据 形式 Data type 格式 Format 数据 来 源 Data sources 
地 形 数据 Topographic data 数字 高 程 ( DEM) .水系 图 GRID 或 .shp file SRTM90M 数据 .全 国 水 系 图 
土地 利用 Land use data 土地 利用 图 GRID 或 .shp file Landsat TM 遥感 影像 解 译 数据 
土壤 数据 Soil data 土壤 图 GRID 或 .shp file FAO 土壤 数据 
EH 量 - 言 与 泪 
气象 数据 Meteorological data 日 降雨 量 . 日 最 高 最 低 气 温 、 晶 .dbf 或 .txt 气象 站 实测 数据 


辐射 量 、 相 对 湿度 、 风 速 数据 
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42°N 


109° 110° 


111°E 


图 5 艾 不 盖 河 流 气 象 水 文 站 点 分 布 图 


Fig.S Distribution of weather stations and hydrologic stations in the Aibugai Basin 


采用 SWAT-CUP 软件 对 模型 参数 的 敏感 性 进行 分 
析 ,选取 CN2 .SOL_AWC .ESCO 等 10 个 参数 进行 调 参 ， 
采用 百灵 庙 水 文 测 站 1986 年 一 1997 年 实测 月 径流 数 
据 进 行 模型 的 率 定 和 验证 ,将 1986 人 1987 两 年 作为 模型 
的 预 热 期 ;1988 年 一 1992 年 为 模型 的 参数 率 定期 ;1993 
年 一 1997 年 为 模型 的 验证 期 。 采 朋 相 对 误差 (RE) 、 相 
关系 数 (R? ) 和 纳什 效率 系数 (BW) 判别 模型 模拟 效果 ， 
模拟 结果 见 表 5、 图 

从 模拟 结果 来 看 ,使 用 SWAT 模型 对 干旱 地 区 内 陆 
型 流域 进行 水 文 模拟 ,模拟 效果 能 够 满足 要 求 , 且 较 好 
地 模拟 年 答 流 变 化 趋势 ,可 被 用 来 进行 水 文 过 程 模 拟 及 
计算 员 利 用 率 定 妤 的 模型 对 人 研究 区 2010 年 土地 利用 条 
伞 忆 的 水 文 过 程 进 行 模拟 ,以 评价 土地 利用 变化 对 生态 
系统 脆弱 性 的 影响 。 


110° 111°E 


图 6 研究 区 子 流域 划分 图 


Fig.6 Division of the sub-basins in research area 


表 5 SWAT 模型 径流 模拟 结果 
TableS The simulation results of SWAT model 
模拟 效果 (多 年 平均 ) Simulation effect( annual average) 


时 期 Period 

RE R? Ens 
率 定 期 (1988 年 一 1992 年 ) Calibration period (1988 一 1992) 一 8.0% 0.76 0.74 
验证 期 (1993 年 一 1997 年 )Validation period ( 1993 一 1997) 6.1% 0.69 0.63 


RE: 相对 误差 Relative error; R?: 相关 系数 Correlation coefficient; Ews: 纳什 效率 系数 Nash coefficient 
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6.00 一 一 实测 值 -一 一 模拟 值 


流量 Flow/(m3/s) 


参考 类 比 全 国 .各 省 区 (市 .县 ) 尤 其 是 西北 干旱 、 半 干旱 地 区 同类 情况 并 


1989 1990 1991 1992 


199 


3 


年 份 Year 


验证 期 


1994 1993 


图 7 SWAT 模型 模拟 值 与 实测 值 对 比 

Fig.7 Contrast of the simulation in SWAT and the measured data 
2.3 艾 布 盖 河 流域 生态 脆弱 性 评价 
采用 本 文 构建 的 干旱 牧区 生态 脆弱 性 评价 体系 对 艾 布 盖 河 流域 1990 年 .2010 年 生态 胸 弱 性 进行 评价 。 


1996 1997 


结合 当地 实际 设置 各 指标 节点 


值 见 表 6。 
表 6 评价 指标 节点 值 
Table 6 The node value of each evaluating index 

指标 名 称 差 (0.2) 较 差 (0.4) 中 等 (066) 较 优 (0.8) 优 (1.0) 指标 方向 
Index Bad( 0.2) Worse( 0.4) Moderate (0.6) Better( 0.8) Good( 1.0) Index direction 
干旱 指数 Drought index 10 8 5 4 3 逆向 
植被 盖 度 Coverage/% 10 20 30 40 50 正 向 

吉安 重 3 
土壤 容重/(g/em ) 1.6 1.5 1.4 1.3 1.2 逆向 
Bulk density of soil 

吉水 革 散 能 
| 0.3 0M5 0.6 0.8 1 正 向 
了 Potential evapotranspiration of soil water 
均匀 度 Evenness index 下 0.8 0.6 0.4 02 道 向 
草 延 度 Contagion index 40 50 60 70 80 正 向 
生态 适宜 性 Ecological suitability 40. 50 60 70 80 正 向 

资源 | 上 党 3 2 
水 资源 利用 模 数 /( 万 m /km ) 名 9 0 道 向 


Water utilization modulus 


由 层次 分 析 法 4AHR) 计算 各 层 权 重 ,突出 人 类 活动 对 生态 环境 的 影响 ,采用 随机 一 致 性 比率 CR 进行 一 
致 性 检验 ,各 权重 矩阵 CR 均等 于 0, 指标 权重 具有 可 接受 的 不 一 致 性 , 表 7。 


目标 层 Tax&et layer 


表 7 评价 指标 最 终 权 重 表 


Table 7 The final weights of evaluating indexes 


准则 层 Criterion layer 


指标 层 mdex layer 


终 权 重 Final weight 


研究 区 生态 脆弱 性 
Eeologicalvulnerability of 


research area 


2.4 评价 结果 


自然 生态 系统 脆弱 性 ”0.4 


人 类 -自然 耦合 系统 脆弱 性 ”0.6 


干旱 指数 0.286 
植被 盖 度 0.286 
土壤 容 可 
土壤 水 燕 散 能 力 0.286 
均匀 度 0.167 
草 延 度 0.167 


E 


EE 0.143 


E 态 适宜 性 0.333 
水 资源 利 月 


日 模 数 0.333 


0.114 
0.114 
0.057 
0.114 
0.100 
0.100 
0.200 
0.200 


研究 区 20 年 间 生态 脆弱 程度 下 降 明 显 , 各 子 流域 生态 脆弱 性 综合 指数 均 出 现 不 同 程度 下 降 ,其 中 9、10、 
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11 12 .15 子 流域 由 1990 年 的 不 脆弱 、 轻 度 脆弱 发 展 为 2010 年 的 中 度 脆弱 ,16 17 子 流域 由 1990 年 的 中 度 脆 
弱 发 展 为 2010 年 的 重度 脆弱 (图 8) 。20 年 间 , 人 研究 区 整体 生态 脆弱 综合 指数 由 1990 年 的 0.70 下 降 到 2010 
年 的 0.59( 图 9) ,生态 脆弱 等 级 由 轻 度 脆弱 变 为 中 度 脆弱 。 


110° 111°B 110° 111°E 
| 不 脆弱 轻 度 脆弱 中 度 脆 弱 。 上 重度 脆弱 。。 极度 脆弱 


图 8 各 子 流域 生态 脆弱 等 级 图 
Fig.8 Ecological vulnerability level of each sub-basin 


3 分析 与 讨论 —e— 1990 —o— 2010 


3.1 研究 区 生态 脆弱 性 分 析 

从 各 子 流域 生态 脆弱 等 级 的 变化 趋势 来 看 ,研究 区 
北部 及 西部 生态 稳定 性 保持 较 好 ,而 生态 脆弱 等 级 下 降 
主要 发 生 在 研究 区 东南 部 。 究 其 原因 ,20 世纪 90 年 代 
初 ,研究 区 内 的 农田 主要 分 布 在 南部 , 随 着 粮食 需求 的 
不 断 加 大 ,农业 区 逐渐 向 北 扩张 8s2000 开始 ,研究 区 逐 
步 实施 禁 牧 、. 轮 牧 政策 ,研究 轿 张 部 的 草地 得 到 了 有 效 
保护 的 同时 ,为 了 支撑 不 断 增 大 的 牲畜 需求 ,包括 灌溉 
饲 草料 地 在 内 的 耕地 面积 进一步 向 北 扩张 ,草地 被 侵占 


及 农田 灌溉 车 确 的 地 人 征 水 过 度 开发 利用 是 研究 区 南部 图 9 研究 区 各 指标 生态 脆弱 性 指数 
生态 脆弱 性 加 剧 的 主要 原 o 从 研究 区 各 评价 指 标 的 Fig.9 Variation of ecological vulnerability index 


变 伐 情况 间 以 看 出 ,干旱 牧区 草地 资源 的 退化 及 减少 造 
成 的 生态 适宜 性 下 降 , 以 及 耕地 面积 扩大 带 来 的 水 资源 过 度 开 发 同时 利用 效率 低下 是 造成 研究 区 整体 脆弱 等 
级 下 陷 的 重要 原因 。 此 外 ,从 景观 格局 来 看 ,研究 区 景观 连通 性 有 所 下 降 ,破碎 化 程度 进一步 加 剧 , 反 映 出 人 
类 活动 对 上 自然 景观 的 不 断 割 裂 。 
3.2 评价 体系 适用 性 分 析 

本 文 针 对 干旱 牧区 生态 脆弱 特点 ,从 气候 土壤、 植被 .景观 格局 .生态 稳定 性 ,水 源 开发 利用 等 方面 人 手 ， 
以 GIS 为 平台 ,建立 了 一 种 基于 SWAT 模型 子 流 域 划分 的 干旱 牧区 生态 脆弱 性 评价 体系 ,对 艾 布 盖 河 流域 
1990 年 .2010 年 生态 脆弱 性 进行 了 评价 。 我 国 西北 干旱 牧区 地 广 人 稀 , 生 态 结构 相对 简单 ,各 类 测 站 与 统计 
资料 支撑 有 限 , 本 评价 体系 紧 紧 围绕 干旱 牧区 “水 - 草 - 畜 ” 关 系 设置 评价 指标 ,将 观测 数据 与 水 文 模型 模拟 相 
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结合 ,统计 数据 与 空间 地 类 相 联 系 , 点 面 结合 ,能 够 较 好 的 反映 出 研究 区 20 年 间 生 态 脆弱 程度 变化 趋势 及 主 
要 问题 。 
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